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Abstract 

We analyze new data for self-avoiding polygons, on the square and 
triangular lattices, enumerated by both perimeter and area, providing 
evidence that the scaling function is the logarithm of an Airy function. 
The results imply universal amplitude combinations for all area moments 
and suggest that rooted self-avoiding polygons may satisfy a g-algebraic 
functional equation. 

1 Introduction 

Self-avoiding polygons (SAP) on a lattice enumerated by perimeter are the 
canonical model of ring polymers, while when enumerated by area as well they 
are widely used as models of vesicles uB, O, O. Despite their widespread 
applicability, n 
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studies |15|, |l2[ |10j, |13j, |lj]. In this letter we give, for the first time, significant 
exact information, subject to a very plausible conjecture. The validity of this 
conjecture is supported by abundant numerical evidence. 

Consider the scaling behaviour of polygon models, whose two-variable gen- 
erating functions are defined by 



applicability, remarkably little exact information is known about their behaviour 
p|, |21|, almost all available information being as a result of numerical 
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G{x 1 q)=Y J P m .,nX m q n , (1) 

m,n 

where p m , n is the number of polygons, equivalent up to a translation, with 
perimeter m, enclosing an area n, on the triangular or square lattice. Sometimes 
horizontal and vertical steps are distinguished, leading to a 3- variable generating 
function G(x, y, q). 

These models, as well as simpler models, display scaling behaviour in the 
vicinity of the multicritical point {(x,q) — > {x~ , l - )) of the form 

G(x,q) ~ G^(x,q) + (1 - q) s F (^p) , (2) 



where F(s) ~ const. /s 7 as (s — > oo) and F(0) ~ const, with 7 = —4. In 
this letter we derive and solve the differential equation satisfied by the scaling 
function F(s). 

In an earlier paper [£2| we discussed the scaling functions appropriate to very 
simple polygon models that satisfy g-linear functional equations, of the form 

N 

G(x, q) = J2 a fc( x ' l) G (<l kx > l) + b (x, q), (3) 

fc=i 

where au{x,q) and b(x,q) are rational functions in x and q. In particular, 
for polygon models such as rectangles, Ferrers diagrams and stacks flgl , the 
functional equation takes the simple form 

G s (x, y, q) = -G s (qx, y, q) + , 

(1 — qx) s 1 — qx 

where s = 0, 1,2 characterizes rectangles, Ferrers diagrams and stacks respec- 
tively. This can be solved by iteration to give 

00 n n 

^'"'^Sd-^d' (4) 

where (t; q) n = IIfc=o (•*• — 1 k t) denotes the g-product. 

For more complicated classes of polygons, such as staircase polygons, column- 
convex polygons, convex polygons, directed and convex polygons, the generating 
functions G(x, q) are known (see [fy, || and references therein). All these solvable 
cases, however, share the property that they are essentially one-dimensional, and 
so functional equations can be constructed that "build up" the polygon from 
smaller pieces. No such construction is possible for SAP. 

For example, staircase polygons, which are defined by the rule that their 
perimeter consists of two directed paths, which proceed only in the positive x 
and y direction, share a common starting point and end point, and otherwise 
do not intersect. They satisfy the following non-linear functional equation |17| : 

G{x,y,q) = - 1 : &{qx,y,q), (5) 

1 — qx 1 — qx 

which allows an iterative solution to be developed. It is 

G{x,y,q) =qxy—- — , (6) 

L,{x,y) 

where L(x,y) is a g-deformed Bessel function 

(-l) n + m x n y m q( n+ ™ +1 ) 



L{x lV )= Y. 



W,q)n{q;q), 



The generating function by perimeter and area, G(x,x,q) was first given in [4|, 
while the generating function by width and area, G(x, l,q) was given in [3|, 
in terms of g-Bessel functions. The full three-variable generating function was 
subsequently obtained H, and in the above x — y symmetric form in pj. In 
[p~7| a semi-continuous version of this model was solved, which has the same 
asymptotic behaviour as the discrete model given here. In pi, the scaling 
function (0) for the discrete model (with <\> = | and 6 — i) has been extracted 
from the solution (g) . 

Below we show how to derive the (non-linear) differential equation satisfied 
by the scaling function for staircase polygons, using the method of dominant 
balance | J22| on the quadratic functional equation (j^). It is 

F(s) 2 ~ 4AF'( S ) - f s = 0. (7) 

with /q = 1/2, /i = —1/32 at the isotropic point y c — x c — 1/4. This is 
a special Riccati equation, the solution being the logarithmic derivative of an 
Airy function. 

Now, earlier numerical work |17| prompted one of us |19| to point out 
that the SAP generating function for rooted self-avoiding polygons G^(x, q) — 
x-j^G(x,q) has the same exponent values <fi = | and 8 = | as do staircase poly- 
gons, and to ask whether the scaling functions might be the same. In this work 
we answer that question in the affirmative. Of course, this is not to suggest 
that rooted SAP satisfy a similar, comparatively simple, quadratic functional 
equation to that of staircase polygons. 

However let us tentatively assume that rooted SAP satisfy an g-algebraic 
functional equation of arbitrary degree N, 

N n 

J2 H ak 1 ,...,k n (x,q)Y[G(q k 'x,q) = b(x,q), (8) 

n=l ki,...,k n i=l 

where 0>ki,...kn{ x iQ) an d b(x,q) are polynomials in x and q. We require the 
sum to be finite. Then, remarkably, by the method of dominant balance, it fol- 
lows that one is inexorably led to the same differential equation for the scaling 
function (R) as found for staircase polygons, though with different constants /o 
and /i. These two constants can be expressed in terms of the amplitude of 
the perimeter generating function G(x, 1) and the amplitude of the first area 
moment A 1 - 1 ' (x, 1) = ^ m np m , n x m of the perimeter generating function. Fur- 
ther, as we show below, the first area moment amplitude can be given in terms of 
the critical point x c , for any lattice! From our extensive enumeration data Jl3[ 
we have already estimated the required critical point and perimeter generating 
function amplitude accurately. 

We can then solve this differential equation, and use it to predict the am- 
plitudes of all area weighted moments, A^ k \x, 1) = Ylm n (n-fcv P" 1 ." 3 ''"' which 
may then be compared to our numerical evaluations, based on extensive enumer- 
ation data, described below. We obtain complete agreement, up to the accuracy 
of our numerical calculations for the first 10 area weighted moments, which gives 



us sufficient confidence to conjecture that the scaling function is exact. Given 
the scaling function for rooted SAP, the scaling function for (unrooted) SAP is 
obtained by integration. 

2 Scaling functions and g-algebraic polygon mod- 
els 

The coefficients in the asymptotic expansion of the scaling function (g) are read- 
ily related to the area weighted moments of the perimeter and area generating 
function (|j) by Taylor's theorem. Expanding G(x, q) about q — 1 gives 

G(x,q) = $> m ,„a:"V = £>(aO(l - q)\ (9) 

m.n k 

where the functions gk(x) are proportional to the area weighted moments, 

gk(x) = '~.|' A^ k '(x, 1). If we make the common assumption that the singu- 
lar behaviour of the coefficients is 

9 { r 9) (*) = , ^ Wk + O ((x c - x)^ +1 ) , (10) 

(x c — X) '* 

where fk is the amplitude, then the assumption of the scaling form (0) constrains 
the exponents to be 7fe = -^-^ (Tj, while the amplitudes fk appear in the 
asymptotic expansion of the scaling function ]22 ] 

n*)=E£- en) 

fc=0 

The exponents 9 and <f> can be extracted from knowledge of the perimeter gen- 
crating function go(x) and the first area weighted moment gi(x). For staircase 
polygons, they can be computed recursively from the functional equation (g), 
giving — | and <f>= f |l^|. 

We next assume that the generating function G(x, q) satisfies the g-algebraic 
functional equation (g), and that the limit q — ► 1 leads to a non-trivial equation 
for the perimeter generating function G(x, 1) with a singularity at i = x c . 
Assuming the existence of a scaling function (g) about x = x c and q = 1, we 
show how its differential equation can be obtained from (g). For brevity, we 
state the argument under simplifying assumptions, which are sufficient for the 
following numerical analysis. The full argument will be presented elsewhere p3| . 

Assume that the critical exponents 6 — h and <f> — | are known. For rooted 
SAP, they can be obtained from numerical analysis of the perimeter generating 
function and the first area weighted moment, as explained below. The method 
of dominant balance consists in expanding the functional equation about x = x c 
and q = 1, thereby introducing the scaling variable s — (x c — x)/(l — q)^. In 
the scaling limit q —> 1, this leads to a differential equation for the scaling 



function F(s). If we assume the scaling form (0), the generating function may 
be expanded about the critical point as 

G{q k x,q) ~ G {re9) {(l-e) k {x c -se^) 1 l-e)+e e F{s)+e 1 -^ e kx c F'{s)+0{e 2 - 2 ' t}+B ) , 

. (12) 
where e = 1 — q. Observing the powers of e that occur in each term, it follows 

that the q-algebraic functional equation (g) is a sum of terms the form 

(e e F) mi (e^+OF')" 12 e m * (se*)™ 4 (13) 

with non-negative integers mi, 777,2,7713, 7774. The differential equation results 
from taking the terms with smallest exponents in e. Terms of order e° vanish, 
while terms of O (e 1 ^ 3 ) give a trivial equation and will therefore vanish, if a 
scaling function exists. Thus, the leading contribution in (H) in the scaling limit 
is of order e 2 ' 3 , which, by power counting results in the differential equation ([7J) 
for the scaling function. An exactly solvable example is staircase polygons (piy. 
The solution of (yft is uniquely determined by the asymptotic behaviour (pTj) 
and given by 

F(8) = -4 A | h Al((A) 2/3 8 ), (14, 

where Ai(#) = i / °° cos(i 3 /3 + tx) dt is the Airy function. The coefficients fk 
in the asymptotic expansion (O) of the scaling function can be written 



fk=c k fifl~ h (A: = 2,3,...), (15) 

where the numbers Ck can be computed recursively from 

fc-i 
2c fc+ i + 6fec fc + ^cic fc+ i_i = (A; = 2,3,...). (16) 

1=1 

We have the initial values c\ — \ and C2 = —5/2. The subsequent values are 
c 3 = 15, c 4 = -1105/8, c 5 = 1695, c 6 = -414125/16, c 7 = 472200, c 8 = 
-1282031525/128, c 9 = 242183775, c 10 = -1683480621875/256. 

3 Generation and analysis of series 

In order to test our predictions, we have generated the full perimeter and 
area generating function for square-lattice SAPs up to perimeter 66 and for 
triangular-lattice SAPs up to perimeter 29. For square-lattice SAPs, we have 
generated the first ten area weighted moments of the perimeter generating func- 
tion up to perimeter 86. The series are available upon request. The method 
used to enumerate square-lattice SAP is an enhancement of the method devised 
by Enting KM in his pioneering work. Recent improvements to the method are 
developed and described in [O, 041. For the triangular lattice, the method and 
program used is precisely that described in |Ic|. 



In it was argued that self-avoiding polygons should admit certain universal 
combinations between the amplitudes of the area weighted moments. Here, we 
give them all in terms of the coefficient amplitudes Ak of the functions gk{%) in 
the expansion (pi). The coefficients of the functions gk(x) grow asymptotically 
as 

[x n ]g k {x)~aA k x- n n^-\ (17) 

where a is a constant such that p m ^ n is nonzero only if m is divisible by a. 
Thus cr = 2 for the square lattice and a = 1 for the triangular lattice M. The 
amplitudes A k are related to the amplitudes fk in (J10|) via 



fk 

xj k v{ lk y 



** = Z^T-T- (18) 



By definition, the absolute values of these amplitudes coincide for rooted SAP 
and for unrooted SAP. The conjecture in [pi reads in our notation that all 
products A k An~ <J k are universal. In particular, from amplitude relations given 
in H and in H we infer A\o = —j^, independent of the lattice! Using this and 
the assumed scaling equation, we obtain from (US) and ([Lq) for even k = 2m 



i 2m- 1 



A AZ?n — L AA 



2m (3m-2)!2 



4m-2 



(6m - 3)! ir m 
_ c 2m (3m -2)! 

47r 3m cr 2m (6m -3)!' 

where the c k are given by ([L6|). For odd k ~ 2m + 1, we obtain 

a Aim _ A2m+1 

Si-2m+lA — C2m+\A 1 



(19) 



(3m)!2 2m 7r r 

C2m+1 



(3 m )\ 7r 3 m +l26ni+2 (7 2 m +l- 



(20) 



We now test these predictions against numerical data obtained for rooted 
SAP on the square and triangular lattices. Note that this tests our entire the- 
ory, as these amplitude relations follow from the scaling function we derived. 
A numerical analysis of the singularity of the area weighted moments confirms 
the behaviour assumed in (Q) with 8 = -| and = | within accuracy given 
by series data. As argued above, we obtain estimates of the coefficients of 
the scaling function from estimates of the amplitudes of the area weighted mo- 
ments. These amplitudes are estimated by a direct fit to the expected asymp- 
totic form. That is to say, we fit the the coefficients to the assumed form 
[x n ]gk(x) ~ x~ n n lk ~ 1 X)j>o a i/ n ^ 1 ^ ' ■ If there is no non-analytic correction term, 
then f(i) = i. A square-root correction term means f(i) = i/2. In all cases, our 
procedure is to assume a particular form for f(i), and observe the extent to 
which it fits the data. With the long series we now have at our disposal, it is 
usually easy to see if the wrong assumption has been made — the sequence of 
amplitude estimates a^ either diverges to infinity or converges to zero. Once the 
correct assumption is made, convergence is usually rapid and obvious. Detailed 



demonstrations of the method can be found in pi, O . In this way, we obtain 
the following results for the leading amplitudes: 



amplitude 


exact value 


square lattice 


triangular lattice 


-AiCT 


.79578XKT 1 


.7957(l)xl0- i 


.796(1) XlO" 1 


A 2 A a 2 


.33595xl0- 2 


.3359(l)xl0~ 2 


.335(l)xl0~ 2 


-A 3 A 2 a 3 


.10025xlCT 3 


.1002(1) xl(T 3 


.100(l)xl0" 3 


A 4 A 3 a 4 


.23755xlCT 5 


.2375(l)xl0- 5 


.237(l)xl0" 5 


-A 5 A*a 5 


.47574xlCT 7 


.4757(1) xl0~ 7 


.475(l)xl0~ 7 


A 6 A 5 a e 


.83663xlCT 9 


.8366(l)xl0- 9 


.83(l)xl0- 9 


-A 7 A e a 7 


.13251xl0- 10 


.1325(l)xl0- 10 


.13(l)xl0~ 10 


A 8 A 7 a 8 


.19242xl0~ 12 


.1924(l)xl0~ 12 


.18(l)xl0~ 12 


-A 9 A*a d 


.25947xl0~ 14 


.2594(l)xl0~ 14 


.25(l)xl0~ 14 


A w Ala w 


.32806xl0" 16 


.3280(l)xl0~ 16 


.31(l)xl0~ 16 



It is clear that the results for the first 10 area weighted moments are in ex- 
cellent agreement with the numerical estimates. On this basis we conjecture 
that the scaling function is correct. Accepting this conjecture, it follows that 
the scaling function F^ r > (s) for rooted SAP is 

F (r) {s) = _£i d ln Ai (JL {2 aA Q f 3 8 ) , (21) 



7TCT ds \ X, 

where we have x c = 0.379052277(1), a = 2, A = 0.2811506(1) for the square 
lattice and x c = 0.240917(1), a = 1, A = 0.263936(1) for the triangular 
lattice. We can check this result further by comparing F^(0) with series 
data evaluated at x = x c , according to (y). We predict from (glj) the val- 
ues F^(0) = 0.394188(1) for the square lattice and F^(0) = 0.476161(1) 
for the triangular lattice, which is in agreement with the numerical estimates 
F( r )(0) = 0.39(1) (square lattice) and F< r >(0) = 0.47(1) (triangular lattice). 
The scaling function F(s) of (unrooted) SAP follows from ( |2l| ) as 

Fts) = — In Ai (— (2<jA ) 2/3 s) (22) 

7TCT \X C J 

with exponents 6—1 and cf> = ^. 



4 Conclusion 

We have given what we believe to be the exact scaling functions for self-avoiding 
polygons in the vicinity of the tricritical point, for both the square and triangular 
lattices. The result follows from the assumption that the two- variable generating 
function for rooted self-avoiding polygons satisfies a g-algebraic functional equa- 
tion. It is possible that even more striking results could flow from exploration of 
this assumption, a topic we are currently investigating. A corresponding result 
for the scaling function in the inflated regime, q > 1 is given in J20[, though the 
SAP structure is largely washed out in this regime. 
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